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Abstract 

We completely solve the problem of enumerating combinatorially inequivalent 
d-dimensional polytopes with d + 3 vertices. A first solution of this problem, by 
Lloyd, was published in 1970. But the obtained counting formula was not correct, 
as pointed out in the new edition of Griinbaum's book. We both correct the mistake 
of Lloyd and propose a more detailed and self-contained solution, relying on similar 
preliminaries but using then a different enumeration method involving automata. 
In addition, we introduce and solve the problem of counting oriented and achiral 
(i.e. stable under reflection) d-polytopes with d + 3 vertices. The complexity of 
computing tables of coefficients of a given size is then analyzed. Finally, we derive 
precise asymptotic formulas for the numbers of d-polytopes, oriented d-polytopes 
and achiral (i-polytopes with d + 3 vertices. This refines a first asymptotic estimate 
given by Perles. 

Acknowledgement. The author would like to thank Giinter M. Ziegler for having posed 
the problem to him and having taken much time to explain thoroughly the combinatorics 
of Gale diagrams and to correct in details a first draft of the paper. 

1 Introduction 

A polytope P is the convex hull of a finite set of point of a vector space R d . If P is not 
contained in any hyperplane of R d , then P is said <i-dimensional, or is called a rf-polytope. 
A vertex (resp. a facet) of P is defined as the intersection of P with an hyperplane H of 
R d such that PC\H has dimension (resp. has dimension d — 1) and one of the two open 
sides of H does not meet P. A vertex v is incident to a facet / if v 6 /. 

This article addresses the problem of counting combinatorially different rf-polytopes 
with d + 3 vertices, meaning that two polytopes are identified if their incidences vertices- 
facets are isomorphic (i.e. the incidences are the same up to relabeling of the vertices). 
Whereas general <i-polytopes are involved objects, ci-polytopes with few vertices are com- 
binatorially tractable. Precisely, each combinatorial ci-polytope with d + 3 vertices gives 
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rise in a bijective way to a configuration of d + 3 points in the plane, placed at the centre 
and at vertices of a regular 2/c-gon, and satisfying two local conditions and a global condi- 
tion. As a consequence, counting combinatorial c?-polytopes with d + 3 vertices boils down 
to the much easier task of counting such configurations of d + 3 points, called reduced 
Gale diagrams. Following this approach, Perles [2, p. 113] gave an explicit formula for 
the number of (combinatorial) simplicial <i-polytopes with d + 3 vertices and Lloyd [3] 
gave a more complicated formula for c(d + 3, d) the number of combinatorial rf-polytopes 
with d + 3 vertices. However, as pointed out in the new edition of Griinbaum's book [2, 
p. 121a], Lloyd's formula does not match with the first values of c(d + 3, d), obtained by 
Perles [2, p. 424]. 

In this article, we both correct the mistake of Lloyd and propose a more complete and 
self-contained solution for this enumeration problem. The following theorem is our main 
result: 

Theorem 1 Let c(d + 3, d) be the number of combinatorially different d-polytopes with 
d + 3 vertices. Then the generating function P(x) = c(d + 3, d)x d+3 has the following 

d 

expression, where <p(.) is the Euler totient function: 

. * 1 xr^ 0(e) , ( 2x 3e \ 1 0(e) , / 1 - x e \ 

P(x) = > In 1 - — + > In 

v ; 1 - x 4e I 1 - 2a; 2e / 1 - x ^ 2e V 1 - 2x e J 

e odd \ \ J / e >! \ / 

x (x 2 - x - 1) (x 4 - x 2 + 1) x (x 8 - 2 x 7 + x 6 + 3 x 3 - x 2 - x + 1) 
2 (1 - xf (2x 6 - 4x 4 + Ax 2 - 1) (1 + xf (1 - xf 

The first terms of the series are P(x) = x 5 + 7x 6 + 31x 7 + 116a; 8 + 379a; 9 + 1133a; 10 + 
3210a; 11 + . . ., i.e. there is one polytope with 5 vertices in the plane (the pentagon), there 
are 7 polytopes with 6 vertices in the 3-D space, etc. 

The mistake of Lloyd, pointed precisely in Section 5, is in the last rational term of P(x). 
Lloyd derived from his expression of P(x) an explicit formula for c(d + 3, d), which does 
not match with the correct values of c(d + 3, d) because of the mistake in the computation 
of P(x). We do not perform such a derivation for two reasons: first, several equivalent 
formulas for c(d + 3, d) can be derived from the expression of P(x), so that the canonical 
form seems to be on the generating function rather than on the coefficients. Second, 
explicit formulas for c(d + 3, d) such as the one of Lloyd involve double summations, 
hence require a quadratic number of arithmetic operations to compute c(d + 3,d). As 
opposed to that and discussed in Section 8, the coefficients c(d + 3, d) can be directly 
extracted iteratively from the expression of P(x) in a very efficient way: a table of the N 
first coefficients can be computed with 0(N\og(N)) operations. Using a mathematical 
software like Maple, a table of several hundreds of coefficients can easily be obtained. 

In Section 7, we introduce the problem of counting oriented rf-polytopes with d + 3 
vertices, meaning that two polytopes are equivalent if they have the same combinatorial 
structure and there exists an orientation-preserving of R d homeomorphism mapping the 
first one to the second one. We establish a bijection between oriented <i-polytopes with 
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d + 3 vertices and so-called oriented reduced Gale diagrams of size d + 3, adapting the 
original bijection so as to take the orientation into account. To our knowledge, this 
oriented version of the bijection was not stated before. The bijection implies that the 
task of counting oriented d + 3- vertex rf-polytopes reduces to the task of counting oriented 
reduced Gale diagrams with respect to the size, which is done in a similar way as the 
enumeration of Gale diagrams. As a corollary, we also enumerate combinatorial d + 3- 
vertex rf-polytopes giving rise to only one oriented polytope. These polytopes, called 
achiral, are also characterized as having a geometric representant fixed by a reflection of 
R d . 

Finally, in Section 9, we give precise asymptotic estimates for the coefficients c(d+3, d), 
c + (d + 3,<f), c~(d+ 3,d) counting (combinatorial) d-polytopes, oriented <i-polytopes and 
achiral <i-polytopes with d + 3 vertices. No asymptotic result was given in Lloyd's paper, 
but Perles [2, p. 114] proved that there exist two constants C\ and c 2 such that Ci-j < 
c(d + 3, d) < C2-j, where 7 is explicit, 7 ~ 2.83. Using analytic combinatorics, we deduce 

from the expression of P(x) that c(d + 3, d) ~ c^r, with c an explicit constant and 7 equal 
to the 7 of Perles, but with a simplified definition. Hence this agrees with Perle's estimate 
and refines it. 

Overview of the proof of Theorem 1. In Section 2.1, we give a sketch of proof of the 
bijection between combinatorial (d + 3)-vertex ci-polytopes and reduced Gale diagrams of 
size d + 3. With this bijection, the enumeration of (d + 3)-vertex rf-polytopes reduces to 
the enumeration of reduced Gale diagrams with respect to the size. 

The scheme of our method of enumeration of reduced Gale diagrams follows, in a more 
detailed way, the same lines as Lloyd. The first observation (see Section 2.2) is that it 
is sufficient to concentrate on the enumeration of reduced Gale diagrams with no label 
at the centre and satisfying the two local conditions (forgetting temporarily the third 
global condition). We introduce a special terminology for these diagrams, calling them 
wheels. As wheels are enumerated up to rotation and up to reflection, they are subject to 
symmetries: Burnside's lemma reduces the task of counting wheels to the task of counting 
so-called rooted wheels (where the presence of a root deletes possible symmetries) and 
rooted symmetric wheels of two types: rotation and reflection, see Section 3. 

After these preliminaries, our treatment for the enumeration of rooted wheels differs 
from that of Lloyd, which relies on an auxiliary theorem of Read, requiring to operate 
in two steps. The method we propose in Section 4.1 is direct and self-contained: we 
associate with a rooted wheel a word on a specific (infinite) alphabet and we show that 
the set of words derived from rooted wheels is recognized by a simple automaton (see 
Figure 3(a)). Under the framework of automata, generating functions appear as a very 
powerful tool providing simple (in general rational) and compact solutions in an automatic 
way. We derive from the automaton an explicit rational expression for the generating 
function of rooted wheels. The enumeration of rooted symmetric wheels is done in a 
similar way, associating words with such rooted wheels and observing that the obtained 
sets of words are recognized by automata. The injection into Burnside's Lemma of the 
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rational expressions for rooted and rooted symmetric wheels yields an explicit expression 
for the generating functions of wheels, given in Section 4.4. Theorem 1 follows after taking 
the global condition (called half-plane condition) into account, which requires only some 
exhaustive treatment of cases, see Section 5. 

2 Gale diagrams of (d + 3)-vertex d-polytopes 
2.1 Gale diagrams 

Following Perles and Lloyd, we define a reduced Gale diagram as a regular 2/c-gon, with 
k > 2, that carries no n- negative labels at its centre and at its vertices, with the following 
properties: 

PI: Two opposite vertices of the 2/c-gon cannot both have label 

P2: Two neighbour vertices of the 2/c-gon cannot both have label 

P3: (half-plane condition) Given any diameter of the 2/c-gon, the sum of the labels of 
vertices belonging to any (open) side of the diameter is at least 2. 

In addition, two reduced Gale diagrams are identified if the first one can be obtained 
from the second one by a rotation or by a rotation and a reflection. The size of a reduced 
Gale diagram is defined as the sum of its labels. The following theorem is essential 
in order to reduce the problem of enumeration of polytopes to the tractable problem 
of counting reduced Gale diagrams. Details of the proof can be found in Griinbaum's 
book [2, Sect. 6.3]. 

Theorem 2 (Perles) The number of combinatorially different d-polytopes with d + 3 ver- 
tices is equal to the number of reduced Gale diagrams of size d + 3. 

Proof (Sketch): Given a rf-polytope P with d + 3 vertices v±, . . . , t>d+3, a matrix Mp is 
associated with P in the following way: Mp has d + 3 rows, the ith row consisting of a 1 
followed by the d-vector-position of the vertex Hence, M P has d + 1 columns, and it 
can be shown that Mp has rank d+1. As a consequence, the vector space V(P) spanned 
by the column vectors (Ci, . . . , C^+i) of Mp has dimension d+ 1, so its orthogonal V(P) ± 
has dimension 2. Let (A±, Aq) be a base of V(P) ± and let A be the (d + 3) x 2 matrix 
whose two columns are (Al,^). Then A is called a Gale diagram of P. The matrix A 
can be seen as a configuration of d + 3 points in the plane, each point corresponding to 
a row of A. The combinatorial structure of P, i.e. the incidences vertices-facets, can be 
recovered from A. However, several Gale diagrams can correspond to the same (combi- 
natorial) polytope. One can perform successive reductions, keeping the same associated 
combinatorial polytope, so that the d + 3 points of the diagram are finally located either 
at the centre or at vertices of a regular 2/c-gon. Giving to the centre and to each vertex of 
the 2/c-gon a label indicating the number of points located at it, one obtains a 2/c-gon with 
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labels characterized by the fact that they satisfy properties PI, P2 and P3. In addition, 
it can be shown that this reduction is maximal, i.e. that the combinatorial types of the 
polytopes associated with two inequivalent (i.e. not equal up to rotation and reflection) 
reduced Gale diagrams are different. □ 

2.2 Gale diagrams and wheels 

A first remark is that properties PI, P2, P3 do not depend on the value of the label 
at the centre of the 2/c-gon. Hence the number g n of reduced Gale diagrams of size n 
is easily deduced from the coefficients counting reduced Gale diagrams of size i with 
label at the centre (such reduced Gale diagrams correspond to so-called non-pyramidal 
polytopes): 

n 

9n = 

i=l 

As a consequence, we concentrate on the enumeration of labelled 2/c-gons (meaning that 
only the 2k vertices of the 2/c-gon carry labels) satisfying properties PI, P2, P3. 

A second remark is that Property P3 is implied by Property P2 if the number of 
diameters is at least 5. As a consequence, we will first put aside Property P3 and focus 
on the enumeration of labelled 2/c-gons satisfying properties PI and P2 (counted up to 
rotation and up to reflections). Such labelled 2/c-gons are called wheels. Wheels with 

2 vertices, even though corresponding to a degenerated polygon, are also counted. The 
enumeration of wheels will be performed in Section 3 and Section 4. By definition of 
wheels, the number of reduced Gale diagrams with no label at the centre is obtained as 
the difference between the number of wheels and the number of wheels not satisfying 
Property P3. The latter term, considered in Section 5, is easy to calculate using some 
exhaustive treatment of cases, because wheels not satisfying Property P3 have at most 4 
diameters. 

3 Method of enumeration of wheels 
3.1 Rooted wheels 

A wheel is rooted by selecting one vertex of the 2/c-gon and by choosing a sense of traversal 
(clockwise or counter-clockwise) of the 2/c-gon. See Figure 1(b) for an example 1 . 

Traversing the 2/c-gon from the selected vertex in the direction indicated by the root, 
one obtains an integer sequence (ai, . . . , a 2 k) satisfying the following conditions: 

SI: For each 1 < i < 2k, and a^ + k) mod 2k are not both 0. 

S2: For each 1 < i < 2k, a« and a(i+i) mo d 2k are not both 0. 

1 On the figures, regular 2fc-gons are represented as 2k vertices regularly distributed on a circle, for 
aesthetic reasons and consistence with the terminology of wheels. 
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(a) A wheel. 



(b) A rooted wheel. 



(c) The associated 
integer sequence. 



Figure 1: Example of wheel and rooted wheel. 



An integer sequence satisfying properties SI and S2 is called a wheel-sequence. The 
size of the wheel-sequence is defined as (ai + . . . + a 2 k). Properties SI and S2 are simply 
the respective transpositions of properties PI and P2 on the integer sequence, so that we 
can identify rooted wheels with size n and k diameters and wheel-sequences of size n and 
length 2k. 

3.2 Burnside's lemma 

Burnside's lemma is a convenient tool to enumerate objects defined modulo the action of 
a group, which means that they are counted modulo symmetries. Let G be a finite group 
acting on a finite set E. Given g e G, we write Fix g for the set of elements of E fixed by 
g. Then the number of orbits of E under the action of G is given by: 



where |.| stands for cardinality. 

3.3 Burnside's lemma applied to wheels 

A wheel with size n and k diameters corresponds to an orbit of rooted wheels with size 
n and k diameters under the action of the dihedral group D 2k . Equivalently, using the 
identification between rooted wheels and wheel-sequences, a wheel with size n and k 
diameters corresponds to an orbit of wheel-sequences of size n and length 2k under the 
action of D 2 k = Z 2 fc x {+, — }, where the action is defined as follows, see Figure 2: 




(1) 



(/,+) • (ai,...,a 2fc ) 
(I, -) • (ai, . . . ,a 2fc ) 



(a i+ i, . . . , a 2 fc, ai, . . . , a/) 
(oi-h, a t , . . . , ai, a 2 k, • • • , o, 2+ i), 
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(a) Rotation action. 



(b) Reflexion action. 



Figure 2: The two cases of action of the dihedral group. 



i.e. (/, +) is a rotation and (I, — ) is a reflection. 

Let us now introduce some definitions. A rotation-wheel is a pair made of a rooted 
wheel and of a rotation of order at least 2 fixing the rooted wheel. Equivalently, it is 
a pair made of a sequence (ai, . . . , a 2k ) and of an element (I, +) with I ^ such that 
(I, +) • (ai, . . . , a 2 fc) = (ai, . . . , a 2 fc). A reflection-wheel is a pair made of a rooted wheel 
and of a reflection fixing the rooted wheel. Equivalently, it is a pair made of a sequence 
(ai, . . . , a 2 fc) and of an element (I, — ) such that (I, — ) ■ (ai, . . . , a 2 fc) = • • • , «2fc)- The 
following proposition ensures that, using Burnside's formula, counting wheels reduces to 
counting rooted wheels, rotation wheels and reflection wheels. 

Proposition 3 Let W n ^ k , R n ,k, Rnk> R n k ^ e respectively the numbers of wheels, rooted 
wheels, rotation-wheels, and reflection-wheels with size n and k diameters. Let W(x,u), 
R(x,u), R + (x,u), and R~ (x,u) be their generating functions. Then 



Proof: As wheels with k diameters are orbits of rooted wheels with k diameters under the 
action of the dihedral group D 2 k (which has cardinality 4k), Burnside's formula yields 



Hence ^ 4kW n<k x n u k = ^ R n , k x n u k + ^ Rn, k x n u k + ^ R~, k x n u k , which yields (2). □ 

4 Enumeration of wheels 



4.1 Enumeration of rooted wheels 

In this section, we explain how to obtain a rational expression for the generating function 
R(x, u) counting rooted wheels with respect to the size and number of diameters. 



4-u— — (x, u) = R(x, u) + R + (x, u) + R (x,u). 



(2) 




( R n,k + Rn,k + R n,k) ■ 
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4.1.1 The word associated to a rooted wheel. 

Let s = (ai, . . . , a 2 k) be a wheel-sequence of size n and length 2k. Associate with s the 
following word: 

'==( *)■('') (*) 

Observe that the length of a is the number of diameters of the associated rooted wheel. 
As each letter of a contains a pair of opposite vertices of the 2/c-gon, the fact that two 
opposite vertices are not both (Property PI or equivalently Property SI) translates into 
the following property: 

a is a word on the alphabet A := N 2 \ 

Now let us detail the translation of Property P2 (or S2) on the word a. First, the alphabet 
A is partitioned into three subalphabets: 




B = <J I M with i > 0, j > ol ,C = j ( M with i > ol ,P= { ( . ) with j > 



• I ' 7 ,y - ^i7— ^17— li • 

kJ/ J IV / J IvJ 

Property S2 is translated as follows: 

• cij and Oj + i are not both for 1 < % < k — 1 <^=^> cr, and o" i+1 are not both in X> for 
1 < i < k - 1 

• a k+i and a^+j+i are not both 0forl<i<A; — 1 ^=^> cr, and o" i+1 are not both in C 
for 1 < % < k - 1 

• dk and Ofe+i are not both <^=^> the pair (a±, a k ) is not in C x V 

• ai and a 2 k are not both ^=^> the pair (ui, cr^) is not in V x C 

Hence a is characterized as a word on the alphabet A that contains no factor CC nor 
factor DP and such that the pair made of its first and last letter in not in C x V nor in 
VxC. 

The size of a letter is defined as the sum of its two integers, and the size of the word 
a is defined as the sum of the sizes of its letters. Hence the size of a rooted wheel is equal 
to the size of its associated word. 

Notice that the generating functions of the three subalphabets £>, C, and V with 
respect to the size are: 
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(a) The generic automaton. (b) Automaton recognizing words 

not containing CC or T>T>. 



Figure 3: Automata associated with words not containing CC or W. 
4.1.2 Generic automaton and its generating functions 

First we explain how to enumerate the words on the alphabet A avoiding the factors 
CC and W. The set of these words is recognized by the automaton represented on 
Figure 3(b), obtained from the automaton of Figure 3(a) by choosing {0} as starting 
state (entering arrow) and {0, 1, 2} as end-states (leaving arrows). We call the automaton 
of Figure 3(a) generic because rooted wheels, rotation- wheels and reflection- wheels will 
give rise to languages on A recognized by slight modifications of this automaton. 

For i G {0, 1, 2} and j G {0, 1, 2}, we denote by dj the set of words accepted by the 
generic automaton that start at state i and end at state j. Let Lij(x, u) be the generating 
function of with respect to the size and length of the word. Looking at the starting 
state and first letter of a word recognized by the generic automaton and ending at 0, one 
gets the following system satisfied by the three generating functions L 00 (x,u), L w (x,u) 
and L 2Q (x, u): 

L o(x,u) = 1 + uB(x)L 00 (x, u) + uC(x)L w (x, u) + uD(x)L 20 (x, u) 
L 10 (x,u) = uB(x)L 00 (x, u) + uD(x)L 2 o(x, u) 
L 2 o(x,u) = uB(x)L 00 (x, u) + uC(x)L 10 (x, u) 

Replacing B(x), C(x) and D(x) by their expressions given in (3), this system becomes 

L 00 (x,u) \ I 1 \ 
L w (x,u) + . 
L 20 (x,u) J \0 J 

Solving this matrix equation, one gets explicit rational expressions for Lq Q (x, u), L w (x, u) 



L 00 (x,u) 
L w (x,u) 
L 20 (x,u) 



(l-x) 2 
ux 2 

i ux^ 

\ WW 



ux ux \ 
l-x l~x \ 



X 

ux 
l-x 



T^x ) 
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Figure 4: Automaton recognizing non-empty words not containing CC or W and not 
ending with C (resp. V) if they start with V (resp. C). 

and L 2 o(x, u), for instance: 

ux 2 (l — x) 



L w (x,u) 



1 — x(3 + u — 3x — ux + x 2 + u 2 x 2 ) 



One can similarly define a matrix-equation satisfied by {L i(x, u), Ln(i, u), L 2 i(x, u)} 
and a matrix-equation satisfied by {L 02 (x, u), Li 2 (x, u), L 22 (x, u)}, from which one gets 
explicit rational expressions for these generating functions. 

4.1.3 Expression of the generating function of rooted wheels 

As we have seen in Section 3.1, rooted wheels with size n and k diameters can be identified 
with non-empty words of size n and length k on the alphabet A, avoiding the factors CC 
and T>T> and such that the pair made of their first and last letter is not in C x T> nor in 
VxC. 

The language of these words is recognized by the automaton represented on Figure 4. 
Hence the generating function R(x, u) counting rooted wheels with respect to the size and 
number of diameters satisfies: 



R(x,u) = uD(x)(L 20 (x,u) + L 22 (x,u)) + uC(x)(L n (x,u) + L w (x,u)) 
+uB(x)(L 00 (x,u) + L 01 (x,u) + L 02 (x,u)) 

Replacing the generating functions on the right hand side by their rational expressions 
yields 

ux (u 2 x 3 — 2 ux 3 + 2 ux 2 — x 3 + 4 x 2 — 5 x + 2) 
(u 2 x 6 — ux 2 — 3 x 1 + x 4 + 3 x + ux — 1) (x — ux — 1) 

4.2 Enumeration of rot at ion- wheels 

As follows from the definition of rotation-wheels and from the identification between 
rooted wheels and wheel-sequences, a rot at ion- wheel corresponds to a pair made of a 
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(a) A rooted wheel of . 



(b) A rooted wheel of 71^ . 



Figure 5: The two kinds of rooted wheels with a rotation-symmetry. 



wheel-sequence s = (ai, . . . , a 2 fc) and of an element I G Z 2 fc\{0} such that the sequence s 
is equal to its /-shift. Writing e for the order of / in Z 2fc (hence e divides 2k), the sequence 
s is characterized by the property that it can be written as e concatenated copies of an 
integer sequence (cti, . . . , ai2k/e)- in addition, it is well-known that for each divisor e of 2k 
there are exactly 0(e) elements / of order e in Z 2 fc. This yields the following lemma: 

Lemma 4 Let be the set of rooted wheels whose wheel- sequence can be written as 
e concatenated copies of an integer- sequence. Let (x, u) be the generating function 
of with respect to the size and number of diameters. Then the generating function 
R + (x,u) of rotation-wheels is: 

R+(x,u) = ^0(e)i? (e) (x,w) (5) 

e>2 

Let e > 2 and consider a rooted wheel of TZ^ e \ so that its associated sequence 
(a±, . . . , a 2 fc) consists of e concatenated copies of an integer sequence a = (a±, . . . , a 2 k/ e )- 
We give a combinatorial characterization of the sequence a by distinguishing two cases: 



The number e of copies is even. In this case, the opposite vertex of at on the 2k- 

gon is ai, see Figure 5(a). As two opposite vertices of a wheel can not both have label 
(Property PI), all integers have to be positive. This condition ensures that two 
neighbour vertices of the 2/c-gon are not both (Property P2). Hence, for r > 1, a rooted 
wheel of with size n and k diameters corresponds to 2r concatenated copies of a 

non-empty sequence of positive integers of size n/(2r) and length k/r. The generating 
function I(x, u) counting non-empty sequences of positive integers with respect to the size 
and length is: 

1 [x, u) = u- 



l-xl-U-r^- l — x(l+u)' 



Hence we obtain: 



u r x 2r 



R {2r \x,u) = I{x 2r ,u r ) = — -. (6) 
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The number e of copies is odd. As e is odd and divides 2k, it also divides k. Hence 
k/e is an integer, that we denote by k'. In this case, for 1 < % < 2k', the opposite vertex 
of CKj on the 2/c-gon is cxu+k') mod 2fc'i see Figure 5(b). In addition, for 1 < % < 2k! , the 
next neighbour of a« on the 2/c-gon is mo d2fc'- As a consequence, the fact that 

(di, . . . , a 2k ) is a wheel-sequence is equivalent to the fact that (a±, . . . , a 2 k>) is a wheel- 
sequence. Thus a wheel-sequence associated with a rooted wheel of 7^ 2r+1 ) corresponds 
to (2r + 1) concatenated copies of a wheel-sequence, so that for r > 1: 

R( 2r+1 \x,u) = R(x 2r+1 ,u 2r+1 ), (7) 

where R(x, u) is the generating function of rooted wheels. 

Finally, equations (5), (6) and (7) yield the following explicit expression of the gener- 
ating function of rotation-wheels: 

R + {x, u) = J2 0(2r + l)R(x 2r+1 , ^ +1 ) + tt 2r \-la + u r) ( 8 ) 

r>l r>l ^ > 

4.3 Enumeration of reflect ion- wheels 

We recall that a reflection-wheel is a pair made of a rooted wheel and of a reflection fixing 
it. It can also be seen as a pair made of a wheel-sequence (a±, . . . , a 2k ) and of an element 
I e Z 2k such that (a x , . . . , a 2k ) = (a>i+i, a t ,. . . ,a 1} a 2k , • • • , a 2+ i) 

Lemma 5 Let R^'^Xju) be the generating function of rooted wheels fixed by the re- 
flection (—1,—) and let R^°~\x,u) be the generating function of rooted wheels fixed by 
the reflection (0,—). Then the generating function R~{x,u) of reflection-wheels is: 

R-(x,u) = u^- (R^-\x,u) + R^-\x,u)) . (9) 
au 

Proof: For k > 1 and I e Z 2 k, we denote by TZ^' k ^ ^ ne se ^ °f ro °ted wheels with size n and k 
diameters whose associated sequence verifies (ai, ... , a 2k ) = (a 1+ i, ai, . . . , a±, a 2 k, ■ ■ ■ , a 2+ i). 
By definition, the set 7Z~ k of reflect ion- wheels with size n and k diameters is given by 

U~ k = uf*o } - Observe that if a wheel sequence is fixed by the action of (/, — ), then 

its r-shift is fixed by the action of (I — 2r, — ). As a consequence, TZ^' k ^ * s m bijection with 

TZn ' k ^ if ^ is even (these cases are those of a reflection fixing two vertices of the 2/c-gon); 

and 7Zn' k ^ * s m bijection with TZ^ k ' ^ if / is odd (these cases are those of a reflection fixing 

no vertex of the 2/c-gon). This directly yields \TZ~ k \ = k {^R^' k ^| + k , from which 

Equation (9) follows. □ 

4.3.1 Enumeration of rooted wheels fixed by the reflection (0, — ) 

Let (ai, . . . , a 2k ) be a wheel-sequence fixed by (0, — ). Then (ai, . . . , a 2k ) can be written 
as (P, ai, ■ ■ ■ , CKfc-i, 7, oi k -\, . . . , «i), see Figure 6(a). Observe that (5 is opposite to 7 and 
that, for 1 < i < k — 1, ctj is opposite to a k -i on the 2/c-gon. Then two cases arise: 
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(a) A rooted wheel fixed (b) A rooted wheel fixed 

by the reflection (0, — ). by the reflection (—1, — ). 



Figure 6: The two kinds of rooted wheels with a reflection-symmetry. 



The number of diameters is odd. In this case, we write r for the integer (k — l)/2. 
Then the fact that two opposite vertices of the rooted wheel do not both have label 
(Property PI) is equivalent to the fact that: 



a :- 



(* 

7/ \®2rJ \a r +i 



is a word on the alphabet A : = N 2 \{(°)}. 

In addition, Property P2 (two neighbours can not both have label 0) translates as 
follows: 

• P and a\ are not both and cej and ccj+i are not both 0forl<i<r — 1 •<=>- a 
contains no factor W. 

• 7 and a 2r are not both and aij and aij+i are not both for r+l<i<2r — 1 -<=>- 
a contains no factor CC. 

• a r and a r+ \ are not both ^=^> the last letter of a is in the alphabet A (already 
implied by Property PI) 

Hence a is characterized as a non-empty word on the alphabet A that contains no 
factor CC nor factor W. The set of words on A satisfying this last property is already 
recognized by the automaton represented on Figure 3(b). However, the first letter of 
a, containing the two fixed points, counts once in the sequence (ai, . . . , a 2 k), whereas 
the other letters count twice. Hence we rather use the automaton of Figure 7, which 
recognizes non-empty words avoiding CC and W, and where the first letter of the word 
is read separately. From this automaton, we get the generating function (x, u) of 

rooted wheels fixed by (0, — ) and with an odd number of diameters: 
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Figure 7: Automaton recognizing words not containing CC or T>T>, and reading the first 
letter separately. 



= uB(x)(L 00 (x 2 ,u 2 ) + L 01 (x 2 ,u 2 ) + L 02 (x 2 ,u 2 )) 
+uC(x)(L w (x 2 ,u 2 ) + L u (x 2 ,u 2 ) + L 12 (x 2 ,u 2 )) 
+uD(x)(L 20 (x 2 ,u 2 ) + L 21 (x 2 ,u 2 ) + L 22 (x 2 ,u 2 )) 
(x + 1) (x 3 + u 2 x 3 — 2 x 2 — x + 2) ux (1 — x 2 + x 2 u) 
(x - 1) (x 6 u 4 - u 2 x 4 - 3 x 4 + x 6 + 3 x 2 + - 1) 



The number of diameters is even. The case of even k is quite similar to the case of 
odd k. We denote by r the integer k/2. Then, the fact that two opposite vertices of the 
2/c-gon do not both have label is equivalent to the fact that 

'").( ai ) (•"•* 

is a non-empty word on the alphabet A := N 2 \{Q}, and that a r (which is self-opposite) 
is non 0. Similarly as for odd k, one can see that the 2/c-gon has no neighbour vertices with 
label iff the word a has no factor CC nor factor T>T>. Hence the conditions for the word a 
(including the fact that the first letter of a is counted once and the other letters twice) are 
the same as for the words considered in the last paragraph, so that the generating function 
of these words is R^d^ ( x ' u ) ■ As a consequence, the generating function Rfven (x, u) of 
rooted wheels fixed by (0, — ) and with an even number of diameters verifies: 

X Jb 

(x + l) 2 (x 3 + u 2 x 3 - 2 x 2 - x + 2) ux 
~x 6 u 4 - u 2 x 4 - 3 x 4 + x 6 + 3 x 2 + u 2 x 2 - 1 

Finally, the relation R(°~\x, u) = R ( ^ d ~\x,u) + R^ v '^(x,u) yields: 

R (Q,-) (x u) _ ^ 3 (* + !) + ^ 3 ~ 2x 2 - x + 2) 

K,) (x-l)(x e u 4 -u 2 x 4 -3x 4 + x e + 3x 2 + u 2 x 2 - 1) 1 ' 
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4.3.2 Enumeration of rooted wheels fixed by the reflection (—1,—) 

A wheel-sequence (ai, . . . , a 2 k) fixed by (—1, — ) can be written as (cti, . . . , . . . , a\). 

For 1 < % < k — 1, the opposite vertex of on the 2/c-gon is atk-i, see Figure 6(b). As in 
Section 4.3.1, two cases arise: 



The number of diameters is odd. In this case, we write r := (k — l)/2. As a r is 
self-opposite and ai is opposite to a 2r +2-i for 1 < i < r — 1 the fact that two opposite 
vertices of the 2/c-gon have not both label is equivalent to the fact that: 

< : f - ) ( 

is a word (possibly empty) on the alphabet A and that a r+1 (which is self-opposite) is 
not 0. It is easily seen that the 2/c-gon has not two neighbour vertices with both label 
iff the word a is empty or starts with a letter in B (because a\ and are neighbour 
to themselves) and contains no factor CC nor factor T>T>. Such words just consist of a 
letter in B followed by a word avoiding factors CC and W. As the set of words avoiding 
factors CC and W is exactly recognized by the automaton of Figure 3(b), we can derive 
the following expression for the generating function R^' \x,u) of rooted wheels fixed 
by (— 1, — ) and with an odd number of diameters: 

RoM^M = (l + u 2 B(x 2 )-(L 00 (x 2 ,u 2 ) + L 01 (x 2 ,u 2 ) + L 02 (x 2 ,u 2 )))-^ 

x 2 u (x — 1) (x + 1) (x 2 + u 2 x 2 — 1) 
x 6 u 4 — u 2 x 4 — 3 x 4 + x 6 + 3 x 2 + u 2 x 2 — 1 



The number of diameters is even. In this case, we write r := k/2. Then the 2/c-gon 
has no opposite vertices both carrying label iff 

a l \ ( a r 

is a word on the alphabet A. It is then easily seen that the 2/c-gon has no neighbour 
vertices both carrying label iff a avoids the factors CC and W and starts with a letter 
in B. As mentioned above, such words consist of a letter in B followed by a word recognized 
by the automaton of Figure 3(b). Hence the generating function Riven~\x,u) of rooted 
wheels fixed by (— 1, — ) and with an even number of diameters is 

R { -^-\x,u) = u 2 B(x 2 ) (L 00 (x 2 ,u 2 ) + L 01 (x 2 ,u 2 ) + L 02 (x 2 ,u 2 )) 

(x 2 — u 2 x 2 — 1) u 2 x A 
M 4 X 6 - u 2 x A - 3 X 4 + X & + 3 X 2 + u 2 x 2 - 1 



Finally, the relation R( \x,u) = R odd ' (x,u) + R eV en (x,u) yields 

/ i n . . ux 2 (ux 4 — u 3 x 4 — ux 2 — 1 + 2 x 2 — x 4 + u 2 x 2 — u 2 x 4 ) 
1 ' J u 4 x 6 -u 2 x 4 -3x 4 + x 6 + 3x 2 + u 2 x 2 -l 



;n) 
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4.4 Expression of the generating function of wheels 

From Burnside's formula and from the expressions of the generating functions of rooted 
wheels, rotation-wheels and reflection-wheels, an explicit expression can be derived for 
the generating function of wheels: 

Proposition 6 Let W(x) be the generating function of wheels with respect to the size. 
Then W(x) has the following expression: 

KJ ^ 4e V (l-2x) 2e J ^ 2e Vl-2a;< 

e odd v v ' ' e>l v 

x (x 2 — x — 1) (x 4 — X 2 + 1) 
+ 2(l-x) (2x 6 -4x 4 + 4a; 2 -l)' ( ^ 

where <f>(.) is Euler totient function. 

Proof: A first easy observation is that W(x) = W(x,u)\ , where W(x,u) is the gener- 
ating function of wheels with respect to the size and number of diameters. Notice also 
that W(x,0) = 0, because a wheel has at least one diameter, so W(x) = ^-(x,u)du. 
Hence the expression of j^(x,u) given in Proposition 3 yields 

W(x) — - J u~ l R(x,u)du + - J u~~ 1 R + (x,u)du + - J u^ 1 R^{x,u)du. 

According to Proposition 5, R~(x,u) = u4- ( j R(" 1 '-)(x,m) + R(°-\x,u)). Hence, the 
explicit expressions of R^°~\x,u) and R^~ 1: ~\x, u) obtained in Section 4.3 yield 

- f u- 1 R-(x,u)du = ](R(- 1 '-\x,1) + R {0 '-\x,1)) 
4 Jo 4 

x (x 2 — x — 1) (x 4 — x 2 + 1) 



2(1 -x) (2a; 6 -4a; 4 + 4x 2 -l)' 
Then, the expression of R + (x,u) given in (8) yields: 

/ u~ 1 {R{x,u) + R + {x,u))du = ^0(2r + l) / u^R (x 2r+1 , u 2r+1 ) du 
Jo r > Jo 

+ £>(2r) /V- U lf T - 2r du. 

^ J 1 — x Zr — u r x 2r 



Writing G(x) := J* u~ 1 R(x,u)du and H(x) := J* u^ 1 1 _"^ a . 2 du, and using the change of 
variable y = u 2r+1 for the rth term of the first sum and y = u r for the rth term of the 
second sum, we obtain 



f u-\R(x, u) + R + (x, u))du = J2 ^ + {x 2r+1 ) + £ ^-H 

^° r>0 r r>l r 
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so that 



\r>0 r>l / 

x (x 2 — x — 1) (x 4 — x 2 + 1) 
+ 2 (1 - a;) (2a; 6 - 4a; 4 + 4a; 2 - 1) 

The integral if (a;) is easy to compute: H(x) = — In (1 — 2 x 2 ) + ln (1 — x 2 ). Using Expres- 
sion (4) of R(x, u) and a mathematical package, one finds G(x) = — In (1 — 4x + 4x 2 — 2 a; 3 ) + 
2 In (1 — x). Then, as observed by Lloyd, ln(l — 4a; + 4a; 2 — 2a; 3 ) = 2 ln(l — 2x) + ln(l — 
2a; 3 /(l — 2a;) 2 ), so that the terms of W(x) can be re-combined into (12). □ 

5 Enumeration of wheels not satisfying P3 

As observed in Section 2.2, wheels not satisfying the half-plane condition P3 have at most 
four diameters. Hence the generating function of these wheels is equal to the difference 
between the generating function of wheels with at most four diameters and the generating 
function of wheels satisfying Property P3 and having at most four diameters. 

The configuration of a wheel is obtained by putting a black disk on each vertex of the 
2A;-gon occupied by a positive label and then by removing the labels. Figure 8 features 
the 13 possible configurations of a wheel with at most 4 diameters. Similarly, Figure 9 
shows the 10 possible configurations of a wheel with at most 4 diameters and satisfying 
P3, where a black disk is rounded if the label of the corresponding vertex must be at 
least 2 in order to satisfy Property P3. For each case on Figure 8 and Figure 9, we can 
calculate the generating functions of rooted wheels, rotation-wheels and reflection-wheels 
corresponding to this configuration, and derive from Burnside's formula the generating 
function of wheels having this configuration. For example, the contribution of the 5th 
case of Figure 8 is ^ (2/(a;) 3 + 4/(a; 3 ) + QI (x 2 ) I(x)), where I(x) = x/(l - x). 

Then the generating function of wheels not satisfying Property P3 is obtained by 
taking the difference between the sum of the 13 contributions of Figure 8 (last column) 
and the sum of the 10 contributions of Figure 9 (last column). Observe that Cases 8, 10, 
12, 13 of Figure 8 exactly match Cases 5, 7, 9, 10 of Figure 9. Hence it is not necessary 
to compute the generating functions of these cases as they disappear in the difference. 
The calculation of the difference yields the following expression of the generating function 
Wps(x) of wheels not satisfying Property P3: 

TIT , , x(x 8 -2x 7 + x 6 + 3x 3 -x 2 -x+ 1) , 1oN 

W P&) = , -TT27, ^5 • ( 13 ) 

(x+1) {1 — x) 

Remark: Lloyd did a mistake in the calculation of Wp^(x). Precisely he forgot to 
subtract the term corresponding to the 8th configuration of Figure 9 in his computation 
of the generating function of wheels not satisfying Property P3 and having 4 diameters 
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Figure 8: The 13 possible configurations of a wheel with at most 4 diameters. A vertex 
has a disk iff its label is positive. Here I(x) := xj (1 — x) stands for the generating function 
of positive integers. 
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Figure 9: The 10 possible configurations of a wheel satisfying the half-plane property and 
having at most 4 diameters. A vertex has no disk if it has label 0, has a black disk if its 
label is positive, and has a rounded black disk if its label must be at least 2 in order to 
satisfy P3. Here I(x) := x/(l — x) stands for the generating function of positive integers 
and J(x) := x 2 /(l — x) stands for the generating function of integers of value at least 2. 
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and 2 vertices with label 0. His presentation also has two typos: the label > 2 at the top 
of the top-right diagram of Fig. 4 (page 129) has to be replaced by > 1 (it seems it is just 
a typo as the corresponding generating function is then correctly calculated). The second 
typo is in the term (2, 3) of page 131, where (X~^ has to be replaced by j^z^i ■ 

6 Proof of Theorem 1 

By definition (see Section 2.1), a reduced Gale diagram with no label at the centre is a 
wheel satisfying the half-plane property P3. Hence, the generating function of reduced 
Gale diagrams with no label at the centre is the difference between the generating function 
of wheels, given in (12), and the generating function of wheels not satisfying Property 
P3, given in (13). Then, as discussed in Section 2.2, the generating function of reduced 
Gale diagrams is the multiplication by 1/(1 — x) of the generating function of reduced 
Gale diagrams with no label at the centre. Finally, according to Theorem 2, the number 
of reduced Gale diagrams of size d + 3 is equal to the number c(d + 3, d) of combinatorial 
ci-polytopes with d + 3 vertices. This yields Theorem 1. 

7 Oriented and achiral d-polytopes with d+3 vertices 

This section deals with the enumeration of oriented ci-polytopes with d + 3 vertices, 
meaning that two polytopes P and P' are identified if they have the same combinatorial 
type and there exists an orientation-preserving homeomorphism mapping P to P' . We also 
introduce here oriented reduced Gale diagrams, meaning that two reduced Gale diagrams 
are identified if they differ by a rotation. 

Theorem 7 Oriented d-polytopes with d + 3 vertices are in bijection with oriented reduced 
Gale diagrams of size d + 3. 

Proof : (Sketch) The sketch of proof is very similar to the one of Theorem 2. Hence we 
keep the same notations, i.e. the matrix M P associated with P and the vector space V(P) 
spanned by the column vectors (Ci, . . . , Cd+i) of Mp. An oriented Gale diagram of P is a 
(d + 3) x 2 matrix whose two columns vectors A\ and A 2 form a base of V(P) ± and verify 
Det(Ci, . . . , Cd+i, Ai, A 2 ) > 0. As mentioned in the proof of Theorem 2, the combinatorial 
structure of P is encoded in A, and also in the reduced form of A. More precisely, if two 
polytopes have the same combinatorial structure and the same orientation, then they have 
the same reduced oriented Gale diagram. In addition, if two polytopes P and P' have 
equivalent (i.e. equal up to rotation only) oriented reduced Gale diagrams, then one can 
deduce from it a continuous deformation of P into P', keeping a polytope equivalent to 
P all the way. In particular, P and P' have the same combinatorial structure and same 
orientation. □ 

This bijection ensures that counting oriented ci-polytopes with d + 3 vertices reduces 
to counting oriented Gale diagrams with respect to the size. This task is done in the same 
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way as the enumeration of Gale diagrams, i.e. we first enumerate oriented wheels and then 
substract oriented wheels not satisfying the half-plane property P3. The only difference 
between wheels and oriented wheels is in the application of Burnside's lemma. Namely, 
oriented wheels with k diameters correspond to orbits of rooted wheels with k diameters 
under the action of the cyclic group Z 2/ t, of cardinality 2k. From Burnside's lemma applied 
to the group Z 2 k, it follows that the generating function W + (x, u) of oriented wheels with 
respect to the size and number of diameters satisfies: 

dW + 

2u— — (x, u) = R(x, u) + R (x, u), 
ou 

where R(x, u) and R + (x, u) are the generating functions of rooted wheels and of rotation- 
wheels. Proceeding in a similar way as in the proof of Proposition 6, the following ex- 
pression is obtained for the generating function W + (x) of oriented wheels with respect to 
the size: 

W+(x) = _V^ln('l- „ 2 t° ) + In (Izf.) 

y J ^ 2e V (l-2x) 2e J ^ e \l-2x e 

e odd v v 77 e>l v 7 

Then, oriented wheels not satisfying the half-plane property P3 are enumerated by 
doing the same exhaustive treatment of cases as in Section 5. For each of the 13 
configurations of Figure 8 and each of the 10 configurations of Figure 9, the associ- 
ated generating function of oriented wheels is calculated using Burnside's Lemma (ori- 
ented formulation). For example the contribution of the second case of Figure 9 is 

g (2J(x) 3 + 4J(x 3 )) = v (i-z) 3 — • The generating function W^(x) of oriented wheels 
not satisfying P3 is the difference between the 13 oriented contributions of Figure 8 and 
the 10 oriented contributions of Figure 9. This yields 

Tir , , . x 11 + 3x 10 -3x 9 -7x 8 + 4x 7 + 4x 6 + 4x 5 + 3x 4 -2x 3 + x 

w U x ) = < —^37, — ^ • 

(x + l) (1 — X) 

Then, the generating function of oriented reduced Gale diagrams is equal to 

^ (w+(x) - W+fc)) . 

see Section 2.2 and Section 6 for an explanation. As oriented Gale diagrams of size d + 3 
are in bijection with oriented <i-polytopes with d + 3 vertices, we obtain the following 
enumerative result: 

Theorem 8 Let c + (d + 3,d) be the number of oriented d-polytopes with d + 3 vertices. 
Then the generating function P + (x) := J2 d c + (d + 3, d)x d+3 has the following expression, 
where 0(.) is Euler totient function: 

l-x \ ^ 2e { {l-2xr) ^ e 

x 11 + 3 x 10 - 3 x 9 - 7 x 8 + 4 x 7 + 4 x 6 + 4 x 5 + 3 x A - 

(x+l) 3 (l-x) 6 
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The first terms of the series are P + (x) = x 5 + 7a; 6 + 38a; 7 + 170a; 8 + 617x 9 + 1979a; 10 + 
5859a; 11 + . . .. 

Observe that a rf-polytope with d + 3 vertices either gives rise to two different oriented 
polytopes or to one oriented polytope. In the first (resp. second) case, the polytope is 
called chiral (resp. achiral). It can be shown that a combinatorial (<i+3)-vertex d-polytope 
is achiral iff one of its geometric representations is invariant under the reflection x\ — > —x\ 
(the proof relies on the fact that achiral d + 3-vertex rf-polytopes are in bijection with 
reduced Gale diagram having a reflection-symmetry). It follows from the definition that 
the generating function of achiral d-polytopes with d+3 vertices is equal to 2P(x) — P + (x) 
where P(x) and P + (x) are respectively the generating function of d-polytopes and of 
oriented ci-polytopes with d+3 vertices. Using the expressions of P(x) and P + (x) obtained 
in Theorem 1 and Theorem 8, this yields the following corollary. 

Corollary 9 Let c~(d+3, d) be the number of combinatorially different achiral d-polytopes 
with d + 3 vertices. Then the generating function P~(x) = ^2 d c~ (d + 3, d)x d+3 is equal to 

(2a; 11 + 4a; 10 - 2 a; 9 - 15 x 8 - 5 x 7 + 23x 6 + 15x 5 - 17x 4 - 14a; 3 + 4x 2 + 5x + l) x 5 
(-1 + xf (2x 6 - 4x 4 + Ax 2 - 1) (x + if 

The first terms are P~(x) = x 5 + 7a; 6 + 24x 8 + 62a; 9 + 141a; 10 + 287a; 11 + . . .. 

8 Complexity of the enumeration 

The complexity model used here is the number of arithmetic operations, where an op- 
eration can be the addition of two integers of 0(N) bits or the division of an integer of 
0{N) bits by a "small" integer, of 0(log(7V)) bits. 

Proposition 10 The N first coefficients counting combinatorial d-polytopes with d + 3 
vertices can be calculated in 0(N\og(N)) operations. 

The N first coefficients counting oriented d-polytopes with d + 3 vertices can be calcu- 
lated in 0(Nlog(N)) operations. 

The N first coefficients counting combinatorial achiral d-polytopes with d + 3 vertices 
can be calculated in O(N) operations. 

Proof: In the proof we concentrate on the complexity of the extraction of the coefficients 
c(d + 3, d) from the expression of P(x) given in Theorem 1 (the cases of c + (d + 3, d) and 
c~(d + 3, d) can be treated similarly). 

Given a generating function f(x), we denote by devwif) the development of f(x) up 
to power x N . To calculate the iV first coefficients c(d + 3, d), it is sufficient to compute 
devN ((1 — x)P(x)), where P(x) = Yld c (d + 3,d)x d+3 . Indeed, the dth coefficient fd of 
(1 — x)P{x) verifies fd+3 = c(d + 3, d) — c(d + 2, d— 1). Hence, once the N first coefficients 
fd are known, it takes O(N) operations to compute the N first coefficients c(d + 3,d), 
that are calculated iteratively using c(d + 3, d) — c(d + 2, d — 1) + fd+3- 
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Multiplying the expression of P(x) given in Theorem 1 by (1 —x) and then taking the 
iVth truncation yields: 

de^((l-x)P(x)) = ]T ^fZ deUjv (G(a; e )) + X;^^ JV (^(a; e )) 

e=l,e o<f<f e=l 
+rfeWAr(K(x)), 

where := - ln(l - 2x 3 /(l - 2x) 2 ), if (x) := - ln(l - 2x) + ln(l - x) and K(x) is an 

explicit rational function. 

As K(x) is rational, its coefficients satisfy a linear recurrence with constant coeffi- 
cients. Hence, finding devN{K(x)) requires 0(N) operations. The development of H(x) 
is explicit, H(x) = J2 n ^(2 n — 1)^™, so that finding devN(H(x)) requires 0{N) operations. 
Let G n be the nth coefficient of G(x). Then xG'(x) = '^ in nG n x n . As xG'(x) is a rational 
function, the coefficients nG n satisfy a linear recurrence with constant coefficients. Hence 
the calculation of devN(G(x)) requires O(N) operations. 

Once devN(G(x)), devN(H(x)) and devN{K{x)) are calculated, it remains to do the ad- 
dition of 4>(e)/(Ae)dev N (G(x e )) for odd e from 1 to N. As dev N (G(x e )) has [N/e\ nonzero 
coefficients, its addition takes 0(N/e) operations. Hence the total cost of the addition is 
J2e=i 0(N/e) = 0(N\n(N)) operations. Similarly the addition of (f)(e)/(2e)dev N (H(x e )) 
for e from 1 to N takes 0(N\n(N)) operations. Finally it costs 0(N\n(N)) operations 
to compute dev n ((1 — x)P(x)). □ 



9 Asymptotic enumeration 

The asymptotic number of combinatorial rf-polytopes, oriented <i-polytopes and achiral d- 
polytopes with d + 3 vertices can be obtained from the explicit formula of their generating 
functions, given respectively in Theorem 1, Theorem 8 and Corollary 9. 

Proposition 11 The numbers c(d + 3, d) and c + (d + 3, d) of combinatorial d-polytopes 
and oriented d-polytopes with d + 3 vertices have the asymptotic form: 

where 7" 1 is the only real root of the equation 1 — Ax + Ax 2 — 2x 3 , 7 ~ 2.839. 

Let a be the unique positive root of 2x 6 — Ax 4 + Ax 2 — 1 andQ(x) := P ^1^^-4x^+2) X " • 
Then the number c~(d + 3, d) of combinatorial achiral d-polytopes with d + 3 vertices has 
the asymptotic form 

c-(d + 3,d) ~ (C+(-l) d C) \ d 
where C := Q(a) w 12.1278, C' := Q(-a) w 0.0346 and A := a' 1 w 1.6850. 

Proof: As for the proof of Proposition 10, we only concentrate on the case of c(d + 3, d), 
(the proofs for c + (d + 3, d) and c~(d + 3, d) can be done with the same tools). We use 
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the framework of singularity analysis to derive an asymptotic estimate of the coefficients 
c(d + 3, d) from the development of P(x) at its dominant singularity (for a generating 
function with nonnegative ceofficients, the dominant singularity is the least real value 
where P(x) ceases to be analytic). From the expression of P(x) given in Theorem 1, 
it is easy to check that the dominant singularity of P(x) is the real value p such that 
1 - 2p 3 /(l - 2p) 2 = 0, i.e. the real root of 1 - Ax + 4x 2 - 2x 3 . In addition, P(z) has the 
following singular development at p, holding in a vicinity {\z — p\ < e} fl {z — p ^ R + }: 



Denote by [x n ]f(x) the nth coefficient of a generating function f(x). As the singular de- 
velopment of P(z) holds in a "Camembert" -neighbourhood of p and as P(z) is aperiodic 2 , 
transfer theorems of analytic combinatorics [1] ensure that 



As c(d + 3, d) = [x d+3 ]P(x), this yields the asymptotic estimate c(d + 3, d) ~ 1T> 
where 7 := p~ l . □ 

The following corollary follows directly from the fact that the growth constant of 
c~(d + 3, d) is smaller than the growth constant of c(d + 3, d): 

Corollary 12 The quantity c~(d + 3, d)/c(d + 3, d), i.e. the probability that a combina- 
torial d-polytopes with d + 3 vertices is achiral, is asymptotically exponentially small. 

For example, c~(d + 3, d)/c(d + 3, d) is less than 10% for d = 10 and less than 0.1% for 
d = 20. 
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